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ABSTRACT 

We investigate the supernova-driven galactic wind of the barred spiral galaxy NGC 7552, 
using both ground-based optical nebular emission lines and far-ultraviolet absorption lines 
measured with the Hubble Space Telescope Cosmic Origins Spectrograph. We detect broad 
300 km s“^) blueshifted (—40 km s“^) optical emission lines associated with the galaxy’s 
kpc-scale star-forming ring. The broad line kinematics and diagnostic line ratios suggest that 
the H a emission comes from clouds of high density gas entrained in a turbulent outflow. 
We compare the H a emission line profile to the UV absorption line profile measured along 
a coincident sight line and find significant differences. The maximum blueshift of the H a- 
emitting gas is ~ 290 km s“^, whereas the UV line profile extends to blueshifts upwards of 
1000 km s“^. The mass outflow rate estimated from the UV is roughly nine times greater 
than that estimated from H a. We argue that the H a emission traces a cluster-scale outflow of 
dense, low velocity gas at the base of the large-scale wind. We suggest that UV absorption line 
measurements are therefore more reliable tracers of warm gas in starburst-driven outflows. 

Key words: galaxies: starburst - galaxies: individual: NGC 7552 - galaxies: evolution - 
ISM: jets and outflows 


1 INTRODUCTION 


Feedback from massive stellar winds and supernova explosions has 
long been identified as a mechanism capable of injecting large 
amounts of energy into the local interstellar medium (ISM) of a 
galaxy, resulting in both the heating and potential removal of the 
gas in the form of a wind ( |Larson|1974) . Supemovae-driven winds 
have myriad profound effects on galaxy evolution, influencing the 
shape of the galaxy lumin osity function ([Benson et al.|2003| l, the 
mass-metallicity relation iFinlator & Dave| 2008^ , and the struc- 
ture of galactic disks ( [Scannapieco et al.|2008^ . Such winds have 
been found to be “ubiquitous” in star-forming galaxies where star 
formation surface densities exceed Y,sfr ~ 0.1 Mq yr~^ kpc“^ 
( |Heckman|2002| l. Under these conditions, rapid star formation re¬ 
sults in a large injection of mechanical energy into the local ISM of 
the galaxy by OB stars, Wolf-Rayet stars, and supernovae. 

Although galactic winds are commonly observed both locally 
and, increasingly, in the high-redshift universe (e.g., |Veilleux et ah] 
1 2005 1 ISteidel et al.||2010^ , it has been difficult to determine how 
much mass these winds remove from their host galaxies. The ex¬ 
tent to which galactic winds affect on-going star formation de- 
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pends highly on the ability of such winds to remove gas from 
their hosts, commonly parameterized as the mass loading factor 
r] = Mout /M*, the rate of mass loss due to an outflow as a fraction 
of the global star formation rate. Measured mass loading factors 
are typically around ?7 ~ 1, but many of these measurements come 
with high uncertainties. |Rupke et al.H2005^ measure mass loading 
factors of p « 0.01 - 1 in ~ 45 starburst-dominated galaxies at 
2 < 0.5. These measurements suffer from order-of-magnitude un¬ 
certainties due to large ionization corrections for Na I absorption 
lines. [Bouche et aT]j2012^ measure rj ^ 2 with uncertainties of 
only a factor of 2 for five galaxies studied via background quasar 
absorption, but warn that it may be incorrect to compare outflows 
measured at large impact parameter to the current level of star for¬ 
mation. At higher redshift, [Pettini et al.| (|200 2^ estimate rj ~ 1 in 
a single Lyman break galaxy at z = 2.73. [Newman et ^ ( |2012bl l 
measure ?7 ~ 2 for high-Esfi? systems at a ~ 2, but these mea¬ 
surements suffer from quoted uncertainties of at least a factor of 3. 
The uncertainties are possibly much larger due to uncertainties in 
the electron density measurement because of low-S/N in the [S II] 
lines. More robust measurements of outflow masses and velocities 
will provide better constraints on mass loading factors. 

Mass loss measurements have typically been easier to per¬ 
form in absorption-line studies, where the column density of the 
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absorbing gas is readily recovered from unsaturated lines. Absorp¬ 
tion studies are well-suited to face-on systems where the galaxy 
serves as a background source. The most common absorption probe 
in the easily-accessible visible wavelength region is the Na I D 
AA5889,5896 doublet, which traces cool, neutral gas (e.g., |Chen| 
|et al.|201^ . But the utility of this line is limited due to its relatively 
low ionization potential (~ 5 eV), leaving it highly susceptible to 
ionization from ultraviolet (UV) sources and making it most useful 
in studies of dense, dusty environments ( |Martin|2005||Rupke et al.| 
|2005| l. Even when neutral Na I is present, ionization corrections 
are highly uncertain, leading to large uncertainties in total gas mass 
calculations. More robust absorption measurements can be made in 
the UV using lines like Mg II AA2796, 2804 or Fe II AA2586, 2600. 
These lines have been used to study galactic winds at moderate red- 
shifts {z ~ 0.5 - 1.5; e.g., |Tremonti et al.|2007[|Weiner et al.|2009^ 
where these lines fall in the visible regime, but in the local universe 
these lines require observations from space telescopes, limiting the 
availability of such data. UV observations carry the additional re¬ 
quirement that observed targets must be UV-bright. 

Emission-line studies have typically been limited to edge-on 
systems where the geometry naturally separates extra-planar emis¬ 
sion from that of the disk. Wind energetics and mass estimates are 
substantially more difficult to measure in edge-on systems, how¬ 
ever, due to the outflow velocities being projected into the plane of 
the sky. Emission line studies in face-on systems are unable to spa¬ 
tially separate nebular lines in disk H II regions from extra-planar 
outflow gas. Wind emission will have distinct kinematic features 
at modest spectral and spatial resolution, allowing for the separa¬ 
tion of the two components in velocity space. Spatially resolving 
the galactic rotation curve minimizes velocity smearing between 
spectral apertures, allowing for more accurate recovery of multiple 
velocity components within a single line profile. Galaxies at low 
inclination are free of strong projection effects, allowing for more 
accurate spatial mapping of an outflow. 

Recent estimates of ionized mass outflow have been made in 
high redshift galaxies (z ~ 2) using rest-frame H a emission lines 
( [Genzel et al.||20lH [Newman et al.|[2bl2a^ . These methods rely 
on decomposing the H a emission line into two velocity compo¬ 
nents, a narrow component and a broad component. This broad 
component is interpreted as being generated by the outflowing 
gas, and can thus be used as a direct measurement of the outflow 
properties. This method has also been used recently in the study 
of high-z AGN outflows ( [Genzel et al.||2014| l and luminous and 
ultra-luminous inf rared galaxies (LIRG , Lm > ULIRG, 

Lir > IO^^Lq; Arribas et al.||20M i. This method is promising 
in these high-redshift studies but requires further testing in the lo¬ 
cal universe, where spatially-resolved observations are possible at 
higher signal-to-noise. 

These H a-based methods for estimating the mass outflow rate 
also require testing against absorption-line methods, as absorption¬ 
line measurements are sensitive to a wider range of column densi¬ 
ties. Absorption-line studies therefore may still be sensitive to gas 
unseen in H a observations. A direct comparison of emission and 
absorption data would provide insight into the relative mass ratios 
of gas traced by these different observational methods. 

In this paper, we analyze the galactic wind in the nearby face- 
on spiral galaxy NGC 7552. We analyze visible long-slit spectra 
to measure the properties and mass outflow rate of the wind using 
nebular emission lines, and investigate the validity of the assump¬ 
tion that the broad velocity component traces the galactic wind. 
We also calculate an independent mass outflow rate using UV ab¬ 
sorption lines, and directly compare these two results. We com¬ 


pare these values to other similar works, including studies of sys¬ 
tems at high redshift. In ® we present a general overview of the 
galaxy NGC 7552. In Owe describe our data sources and de¬ 
tail our analysis methods. We examine the results of our analy¬ 
sis in and and discuss broader implications in Finally, 
we summarize our conclusions in 0 For this work we adopt 
Hq = 73 km s”'^ Mpc~^, fim = 0.27, and Ua = 0.73. 


2 NGC 7552 


NGC 7552 is a nearby barred spiral galaxy in the Grus Quartet, a 
group of four large spiral galaxies in the southern sky. It lies at an 
estimated distance of 22.5±1.6 Mpc based on a three-attractor flow 
model (1" = 0.109 kpc; Mould et al.[2000 1 , with a heliocentric ra¬ 
dial velocity of 1586 km s“ ( Pan et al.|2013 i. The galaxy is seen 
largely face-on at an inclination of only 23 - 28° ([Feinstein et al.[ 
[T^ [Pan et al.[[20T^ . Although the morphology of NGC 7552 
does not appear highly disturbed, H I observations of the galaxy 
group show numerous tidal features indicative of interactions be¬ 
tween the group members ( [Freeland et al.|200^ . NGC 7552 has a 
galactic stellar mass of logM* = 10.52, measured by the Spitzer 
Survey of the Stellar Structure in Galaxies (S^G) via a combi¬ 
nation of the 3.6 ^m and 4.5 pm imaging bands on the Spitzer 
Space Telescope jSheth et al.|[2010] (. [Moustakas et al.[ P010[ ( re¬ 
port a metallicity of 12 -|- log(0/H) = 9.16 ± 0.01 using the 
“strong line” photoionization scaling, suggesting the galaxy is sub¬ 
stantially super-solar (Z ~ 3Zq). They also report measurements 
tied to the empirical metallicity scale that are 0.81 dex lower, de¬ 
cidedly sub-solar. However, we have found in our analysis that a 
twice-solar-metallicity stellar continuum model is a good match to 
our data (see § |3.1[[T2l l. Despite the ambiguity in measurements 
of 12 + log(0/H), we believe that the galaxy is of at least solar 
metallicity. 

NGC 7552 has been classified as a “low-ionization nuclear 
emission-line region” (LINER) based on its [O I] emission ([ Durret[ 
[& Bergeron] I988|l, raising questions of whether it hosts an active 
galactic nucleus (AGN). [Grier et ah] ( [2011^ measure a nuclear X- 
ray luminosity of Lo.a-s kev = 1-73 x 10^'^ erg s“^, corrected 
for our assumed distance. This is well-below the X-ray luminosity 
of nearby low-luminosity AGN (e.g., [Asmus et al.|2011[ l and thus 
we infer that if an AGN is present, it is a minor contributor to the 
galaxy’s bolometric luminosity. 

The galaxy just barely qualifies as a LIRG with a tot al infrared 
luminosity of Lm = 1.2 x IO^^Lq (Sanders et al. 2003 again cor¬ 


rected for our assumed distance). Using the calibration of [Calzetti| 
( [2013^ , which assumes a Kroupa IMF ( [Kroupa[|200I[ l, this lumi¬ 
nosity corresponds to a star formation rate of M* = 13 Mq yr“^. 
Another approach is to use both the total infrared luminosity and 
the far-ultraviolet luminosity in order to account for both ob¬ 
scured and unobscured star formation. Using measurements from 
the Galaxy Evolution Explorer (GALEX; [Martin et al.[[2005|l, we 


adopt a far-UV luminosity of Lpuv = 2.38 x IO^Lq (Gil de 


[Paz et aL| [2007[(. Combined with Lm from above, and using the 
IR-(-FUV calibrations of [Hao et al.[ ( [2011) , this measurement yields 
M* = 9.9 Mq yr“^. Va rious other indica tors yield global SFRs 
from 15-47 Mq yr“^ jPan et aL]|2013 and further references 
therein). 

One of the defining morphological features of NGC 7552 is 
its circum-nuclear starburst ring (e.g., [Feinstein et al.|1990|[Forbes[ 
[et aL|1994^ . This ring is approximately 0.5 kpc in radius and con¬ 
sists of numerous “knots” of heavy, embedded star formation. The 
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Rest-frame wavelength (A) 

Figure 2. An example of the stellar continuum fit. A representative spec¬ 
trum is shown in black and our stellar continuum model is shown in orange. 
The three blue shaded regions denote the spectral windows that were used 
to match the model to the observed spectrum. 


Figure 1. An i?-band image of NGC 7552 from the SINGS project overlaid 
with our eight O'.'8-wide slit positions in orange. The white dashed rectangle 
denotes the 4 kpc X 3 kpc region detailed in Figurej^ 


SFR within this ring alone is as high as 10-15 Mp yr with a 


total infrared luminosity of 5.5 x 10 ^°I/q jPan et al. 


2013 


and fur¬ 
ther references therein). Using near- and mid-infrared observations, 
[Brandi et al.| ( |2012l > identify nine prominent star clusters within the 
ring, with stellar ages ranging between 5.6 Myr and 6.3 Myr. 


3 DATA 
3.1 Visible 

Our main data set consists of 16 long-slit spectra obtained at eight 
different position angles, each centered on the nucleus, using a 
0”8-wide slit. The eight position angles range from 35°-185°, 
providing coarse two-dimensional spatial information across the 
face of the galaxy (see Figure [TJ. Our data were obtained us¬ 
ing the 100-inch Irenee du Pont telescope at Las Campanas Ob¬ 
servatory using the Modular Spectrograph instrument. The CRAF 
CCD with 12 /rm pixels was used, giving a native plate scale of 
0”292 pix~^ before applying 2x2 on-chip binning. The observa¬ 
tion dates span 1996 August 15-18. Our observed spectral band¬ 
pass is 6400-6900 A, covering the rest-frame [N Ilj/Fl a complex 
and the [S It] AA6716,6731 doublet. Sky line measurements indi¬ 
cate a velocity resolution of approximately 40 km s“^ full-width- 
at-half-maximum (FWFIM). 

The data were reduced using standard IRAF long-slit spectral 
reduction techniques. The data do not include separate sky point¬ 
ings, so we created sky apertures for each observation by identi¬ 
fying regions that were free from nebular emission under visual 
inspection. The median of the rectified sky spectra was subtracted 
from the data. We note that in general there are no strong sky lines 
in our bandpass. We then removed these sky apertures from our 
subsequent analysis. We do not apply a correction for telluric ab¬ 
sorption, as the prominent telluric features do not overlap our lines 
of interest. 


In order to remove the stellar continuum from our observa¬ 
tions, we subtracted a model starburst continuum from each aper¬ 
ture. We find that a twice-solar-metallicity, 100-Myr constant-SFR 
model from jGonzalez Delgado et J^ ( |2005| ) adequately matches the 
continuum over our target spectral range. This model accounts for 
deep stellar H a absorption, providing an important correction for 
lines with low equivalent widths (EW < 10 A). The model has a 
wavelength spacing of 0.3 A (~ 13 km s“^), and we then con¬ 
volve the model by our instrumental resolution and approximate 
narrow-component line width (~ 100 km s“^; see to account 
for stellar velocity dispersion. To match the model continuum to our 
data, we first used the median values within three different spectral 
windows to create second-order polynomial representations of both 
the model and the data. We then scaled the model to the data using 
the ratio of the two polynomials as a pixel-to-pixel scaling factor. 
This scaled model was then subtracted from the data to remove the 
continuum flux. A representative example fit is shown in Figure]^ 
Prior to fitting the emission lines, we adaptively binned the 
spatial apertures to create apertures with S/N(H a) > 3. This only 
proved necessary outside the inner 1 kpc where the galaxy’s ro¬ 
tation curve has become fairly flat, greatly reducing the amount 
of velocity smearing introduced by binning these apertures. We 
performed line fitting using a custom wrapper to the MPFIT pack¬ 
age ( |Markwardt|2009^ , a non-linear least squares fitting code writ¬ 
ten in the Interactive Data Language (IDL). This custom wrapper 
simultaneously fits Gaussians to the five nebular emission lines 
within our spectral window: [N II] AA6548,6583, H a, and [S II] 
AA6716,6731. We ultimately fit the data in two passes. We initially 
fit our data with a single Gaussian component for each line. A sin¬ 
gle line component did not produce a good fit to the data, moti¬ 
vating the addition of a second Gaussian component to the fit. For 
the second pass we simultaneously fit two Gaussian components to 
each emission line. We used the results of the single-component fit 
as a “first guess,” providing initial parameters for approximate line 
centers, widths, and fluxes. For our fitting we restrict the line widths 
of the narrow and broad components (Gaussian a < 70 km s“^ and 
71 < a < 300 km s“^, respectively). We experimented with un¬ 
restricted line widths in our fitting and found that this only has a 
~ 4% effect on our calculations in sQ as such, we retain the line 
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Figure 3. Line profile fits for a representative spectral aperture of our data. The observed spectrum is in black, the narrow and broad components are shown 
in dashed-blue and orange, respectively, while the purple dot-dashed line shows the total tit, or the sum of the two components. The dashed gray lines in the 
residual plots denote the ztlcr levels. Left: The [N II]/H a complex. Right: The [S II] 6716/6731 doublet. 


width restrictions to keep the narrow component and broad compo¬ 
nent easily-separated in our analysis. 

For each velocity component, the flux of each spectral line 
is allowed to vary, while we assume each spectral line shares the 
same velocity center and line width. We do, however, fix the [N It] 
flux ratio at/essa^esrs = 3 in accordance with atomic physics 
( |Osterbro^|1989t . Line widths are also corrected for instrumen¬ 
tal resolution (cjinst = 17 km s“^). It is possible that the various 
emission lines do not share the same kinematics. However, we gen¬ 
erally find a good two-component fit (Xred ~ 1) in a large ma¬ 
jority of apertures, and find that our data do not support additional 
velocity components. Our fitting is also as-restrictive or often less- 
restrictive than many similar recent studies (e.g., [Newman et al.| 
|2012b||G^el et al.|2014[|Ho et al.|20l'4t . 

An example of the result of our line fitting is shown in Fig- 
ure[^ Our fitting routine produces good fits in almost all measured 
spectral apertures. A small number of apertures within ~200 pc of 
the galactic nucleus show larger fit residuals. In apertures where 
Xred > 2 and S/N(H a) > 10, we experimented with adding an 
additional broad Gaussian component to the fit. However, this ad¬ 
ditional component in the fit proved to be unstable and did not give 
consistent results. We therefore conduct our analysis solely on the 
basis of a two-component fit. 

Because our analysis is concerned with spatially-resolved 
properties of the galaxy, we must ensure that our observations are 
properly registered with respect to each other. Our data files do 
not have WCS information recorded in the headers, so a manual 
alignment was required. We have no spatial information perpen¬ 
dicular to the slit direction, so we only investigate the centering of 
the galaxy along the length of the slit and assume the slit passes 
through the nucleus. We made an initial estimate of the galaxy cen¬ 
ter location in each slit by selecting the aperture with the highest 
continuum flux. In order to determine any required offsets from 
our initial estimate, we compare the continuum levels of our spec¬ 
tra to radial profiles extracted at each of our slit position angles 
from archival /?-band imaging from the SINGS project (|Kennicutt| 
|et al.|200^ . This comparison showed good slit alignment for most 
slit positions, and only small corrections (1-3 pixels) were needed 
for two slit positions. As a check of our registration, we compared 


all the median continuum values measured at r = 0 and found a 
standard deviation of 6 percent across all slits. The largest devia¬ 
tion from the mean of the median continuum levels is 10 percent. 
Some slit positions show evidence of requiring sub-pixel shifts for 
ideal registration, but we have chosen not to interpolate our data on 
the sub-pixel level as our analysis does not strongly depend on the 
precise 2D registration of our observations. 

To further check our registration, we created an interpolated 
2D map of our measured total H a emission and compared the re¬ 
sult to an H a narrowband image, again from the SINGS project. 
We used the SurG^ interpolation software package for our 2D 
image reconstruction. The comparison is shown in Figure]^ The 
reconstruction shows good qualitative agreement with the narrow- 
band image data. In particular, we successfully recover the well- 
known star-forming ring of NGC 7552. 

One goal of this work is to calculate the mass outflow rate of 
NGC 7552 in a manner consistent with recent studies performed 
at high redshift {z ~ 2; e.g., [Newman et al.||^12b) . This cal¬ 
culation depends on the observed H a broad-line luminosity. To 
ensure accurate absolute spectrophotometric calibration, we have 
matched our integrated H a counts to calibrated H a narrow-band 
imagery from the SINGS project. This image was acquired through 
a narrow-band H a filter with a central wavelength of 6596 A 
(6561 A in the galaxy’s rest-frame) and an 18 A FWHM, result¬ 
ing in negligible [N II] contamination. The broadband image has a 
calibration accuracy of approximately 10 percent. 

3.2 Ultraviolet 

We also use an archival ultraviolet (UV) spectrum of NGC 7552 
obtained with the Cosmic Origins Spectrograph (COS) instrument 
( [Osterman et al.[[2011[ [Green et al.[[2012) on-board the Hubble 
Space Telescope (HST). The spectrum was acquired as part of HST 
program 12173 (PI: Leitherer), which explored the UV properties 
of four infrared luminous galaxies. [Leitherer et al.|p013^ analyzed 
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Figure 4. H o narrowband image of NGC 7552 from the SINGS project 
(top) and our 2D reconstruction of our total observed H a flux (bottom). 
Both images use logarithmic color scaling. The images show the central 
4 kpc X 3 kpc region of the galaxy highlighted in Fig. [T] The black circle 
represents the COS primary science aperture (2'.'5 diameter) at the pointing 
of our archival UV absorption data. For reference, we also show a longslit at 
PA = 125°. Data from this slit position are shown in Fig.[^ The crosshairs 
mark the nucleus of the galaxy. 


the stellar content, Lyman-o? profiles, and interstellar medium ab¬ 
sorption lines of all four galaxies and estimated a mean mass out¬ 
flow rate for the sample. Here we take a closer look at the multi¬ 
phase gas in NGC 7552 and focus on the comparison of outflow 
properties inferred from the UV and H a lines. 

The COS spech'um was acquired using the 2''5-diameter Pri¬ 
mary Science Aperture with the G130M grating. Fortuitously, the 
COS pointing falls along one of our long-slit positions, allowing 
for a direct comparison of the same physical location in the galaxy 
(see Figure]^. The data cover the wavelength region 1139.52 - 
1451.53 A. The point source spectral resolution of COS/G130M 
is i? « 17,000 (18 km s“^). Measured line widths of unre¬ 
solved Milky Way foreground lines in the spectrum indicate that 
the achieved resolution is closer to 100 km s^'^FWHM. The lower 
resolution is a consequence of the fact that COS is optimized for 
point sources and NGC 7552 fills the aperture. Additional data ex¬ 
ist from the Faint Object Spectrograph (FOS) and the Space Tele¬ 
scope Imaging Spectrograph (STIS) instruments, both also aboard 
HST. A visual inspection of this additional data shows qualitatively 


similar features to the COS spectrum, but at substantially lower res¬ 
olution; we therefore use only the COS data in our analysis. 

In order to minimize the impact of micro-channel plate de¬ 
tector fixed-pattern noise, NGC 7552 was observed at two cen¬ 
tral wavelength positions shifted by 18 A for a total of 5599 sec¬ 
onds ( [Leitherer et al.|[2013^ . The individual spectra were down¬ 
loaded from the MAST server and processed through the CalCOS 
pipeline, version 2.20.1. The individual exposures were co-added 
using a procedure developed by |Wakker et al.| ( [2015| l to optimally 
align the spectra. Strong absorption lines in individual exposures 
are cross-correlated to give a wavelength-dependent wavelength so¬ 
lution, and then the raw counts are combined and converted into 
flux values. The errors are calculated from the raw gross counts, 
after accounting for the proper background levels (as defined in 
the COS handbook ( [Holland et al.|2014l l). Special attention is paid 
to edge effects and fixed pattern noise features, which are known 
to affect the COS data ( jPanforth et al.|[2010[ . We correct the fi¬ 
nal spectrum for foreground Milky Way dust attenuation using 
E{B — V) = 0.012 i Schlafly & Finkbeiner||201 1 1 and a Milky 
Way dust attenuation curve from ( |Fitzpatrick|1999^ . 

Our goal is to use the rich array of interstellar medium (ISM) 
absorption lines in the COS spectrum to probe gas in the outflow. 
To accurately measure the ISM lines, we must first model and re¬ 
move any features in the background light source: the stellar contin¬ 
uum. We use the standard technique of UV spectral synthesis (e.g., 
[Tremonti et al.|200T) to identify a best fit model of the stellar con¬ 
tinuum. We use stellar population models produced by Starburst99 
(SB99; [Leitherer et al.|1999[ with a Kroupa IMF, a power-law in¬ 
dex of 1.3 for the low mass slope, an index of 2.3 for the high 
mass slope, a high mass cut-off at 100 Mq, and the Geneva stellar 
evolution tracks for non-rotating stars with high mass-loss. We uti¬ 
lized the high-resolution, fully theoretical model atmospheres cal¬ 
culated using the WM-Basic code ( [Leitherer et al.|20l0| l instead of 
the default empirical libraries. The theoretical libraries provide an 
excellent match to the spectral resolution of our data and are free 
of the ISM absorption lines that contaminate the empirical stellar 
libraries. 

We took an iterative approach to precisely determine the red- 
shift of the COS spectrum. This was necessary because the spec¬ 
trum is of an off-nuclear position on a steeply rising portion of the 
rotation curve. We used the galaxy redshift as a starting point, found 
the best-fit SB99 continuum model, and then cross-correlated the 
model with the data (with ISM lines masked out) to find an im¬ 
proved redshift. The data was then re-fit with the SB99 models and 
the procedure repeated. 

We fit the COS spectrum with a linear combination of SB99 
single-age stellar population (SSP) models with ages ranging from 
1-20 Myr. We tested models with different metallicities and ulti¬ 
mately adopted models with a twice-solar abundance because they 
provide a better fit to the Si IV stellar wind lines. Our best-fit model 
is dominated by a 3 Myr-old stellar population with small contri¬ 
butions from 1, 4, and 15 Myr-old stars. The mean light-weighed 
age is 3.4 Myr, slightly younger than the ages measured for heavily 
embedded star clusters in the nuclear ring (5-6 Myr;|Brandl et al.j 

I20T2I 1. 


Our stellar population model fits include the effects of dust us¬ 
ing a [Calzetti et al.[ ( |2000[ ( dust attenuation law. The best-fit stellar 
continuum reddening is j5(73 —U) = 0.64, approximately 0.1 mag 
higher than the value found by [Leitherer et al.| ( [2013[ l using simi¬ 
lar methods. This small discrepancy in the reddening is indicative 
of the level of systematic uncertainty. It largely reflects differences 
in the assumed star formation history: [Leitherer et al.[used a SB99 
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model with a continuous star formation history, whereas we use a 
sum of SSP models. 

Our final best-fit dust-reddened SB99 model is compared to 
the data in Figure The observed spectrum is divided by the 
model to produce a normalized spectrum which is used for ana¬ 
lyzing NGC 7552’s ISM lines (see A final step in our data 
processing was to fit and remove the ISM absorption arising from 
the Milky Way. This was carried out using the procedure described 
in |Chisholm et al.|f2014[ >. Fortunately, the Milky Way ISM lines are 
blueshifted by approximately 1550 km s“^ relative to the systemic 
velocity, and thus for the most part do not overlap with absorption 
associated with NGC 7552’s ISM. 


|20I I| >, while we see consistent broad line widths across a region 
~ 2 kpc in diameter. AGN-driven winds have been invoked to ex¬ 
plain spatially-extended broad H a emission lines in the central 
regions of some z ~ 2 galaxies ( [Genzel et al.|2014l l. However, we 
consider this explanation unlikely considering the low luminosity 
of NGC 7552’s nuclear source. The tight correspondence between 
the broad lines and the inner star-forming ring suggests that the 
broadening mechanism is related to the star formation process. It 
is unlikely that these broad line widths result from typical galactic 
H II regions, and we therefore consider a number of other possible 
origins. 


4 MEASURING THE OUTFLOW USING OPTICAL 
EMISSION LINES 

Our fitting routine successfully identifies two emission line com¬ 
ponents over much of the galaxy, with the most significant detec¬ 
tion of broad emission occurring within 1 kpc of the galaxy cen¬ 
ter. Representative examples of our extracted emission line data are 
shown in Figure]^ showing data from PA = 125°. In Figs.|^a) 
and (b) we see a clear separation between the narrow and broad 
components in both the velocity offsets and line widths of the com¬ 
ponents. The line centers of the broad components are consistently 
blueshifted by An « 30 — 50 km s “ ^ compared to the narrow com¬ 
ponents, but still largely trace the rotation curve seen in the narrow 
lines. The line widths of the broad components also show a distinct 
separation from the narrow components; we consistently measure 
line widths of ~300 km s“^ FWHM for the broad components, 
compared to ~100 km s”^ FWHM for the narrow components. 
In Fig. I^c) we show the star formation surface density T^sfr 
for both velocity components, calculated from the H a luminosity 
and uncorrected for extinction, using the relation SFR (Mp) yr~^) 

j2m 


Hao et al. 


. The hori- 


= 5.53 X (erg s“^) from 

zontal dotted line denotes the threshold for driving outflows from 
|Heckman| \2002\ . Fig. |^d) shows the ratio of broad H a flux to 
narrow-line H a flux, /b//n- The broad H a flux is consistently 
at least half that of the narrow flux, and the ratio reaches values 
above one in regions of highest total H a luminosity. Figs. |^e) 
and (f) show the emission line ratio diagnostics log([S II]/H a) and 
log([N II]/H a). Overall, we observe larger line ratios in the broad 
emission line components in both diagnostics with depressed line 
ratios in regions of highest total H a luminosity. 


4.1 Line broadening mechanisms 

As Figure]^ shows, the line widths of the narrow components are 
around 100 km s“^ FWHM (cr ~ 40 km s”'^), broad enough to be 
resolved by our observations. These line widths are consistent with 
those seen in H II regions of other LIRGs (e.g., |Rich et al.|2011 [ 
|Arribas et al.|2014} , and a factor of ~ 2 larger than generally seen 
in local spiral galaxies jEpinat et al.|2010y The spatially-resolved 
velocity centers of the narrow lines show clear signs of galaxy rota¬ 
tion, further supporting the interpretation that these lines arise from 
star-forming H II regions embedded in the disk of the galaxy. The 
broad components have widths of ~ 240 — 350 km s“^ (a « 100 
- 150 km s“^), substantially broader than the narrow components. 

Even if NGC 7552 hosts a low-level AGN, the broad line 
widths are not a result of AGN activity. Reverberation mapping 
has shown the diameter of the broad line regions (BLR) of AGN 
to be much smaller than a parsec ^entz et al.|2009| [Brewer et al.| 


4.1.1 A starburst-driven outflow 


A number of recent works (e.g., Genzel et al.|201 l[|Soto et al.|2012t 


[Newman et al.|2012a[ [Arribas et al.|2014 1 suggest that broad op¬ 

tical emission lines are associated with large-scale galactic winds. 
This conclusion is based in part on the finding of predominantly 
blueshifted line profiles. And as Eigure|^a) shows, in NGC 7552 
the broad line velocity centers are consistently blueshifted by 30 
- 50 km s“^ with respect to the narrow emission across the entire 
star-forming ring (|r| < 0.5 kpc). The ring shows very high lev¬ 
els of extinction with an average Ay = 5.3 prandl et al.|201^ , 
implying that H a emission originating on the far side of the disk 
would be attenuated by a factor of ~ 50. Thus the blueshifted broad 
lines in NGC 7552 cannot arise from in-falling material on the far 
side of the disk; they must be associated with an outflow. 

Outflows are frequently modeled as spherical shells of mate¬ 
rial, moving at a constant velocity Vah. In this simple model, pro¬ 
jection effects lead to observed velocities that are smoothly dis¬ 
tributed between —Vah and +Vsh, with a mean velocity of zero 
and a FHWM of 2vsh. This simple model can be adapted to a bi- 
conical geometry by restricting the angular covering factor of the 
shell. When viewed along the axis of the bi-cone, this leads to re¬ 
duced emission at zero velocity, but the same mean velocity (zero) 
and FWHM (2vsh)- The observed velocities may be affected by ex¬ 
tinction associated with the galaxy disk or with the shell itself. This 
can lead to a reduction of the flux at positive velocities and a net 
blueshift of the line profile. In this physical picture, the observed 
line width is a consequence of projection effects and extinction, 
and the most blueshifted gas traces the true deprojected velocity of 
the shell. 

As we show in Figure |^a), the broad components of 
NGC 7552 have FWHM ~ 240 - 350 km s"^ and blueshifts of 
30 - 50 km s“^ with respect to the narrow components across the 
entirety of the star-forming ring. The simple shell model with ex¬ 
tinction has difficulty reproducing these numbers. If heavy disk 
extinction removes the redshifted component of the line profile, 
then the line width is FWHM ~ Vah and the line centroid is 
Veen ~ Vshl2. Thus, the line width is roughly twice as large as 
the line centroid velocity: FWHM/Ucen ~ 2. In contrast, we mea¬ 
sure FWHM/wcen ~ 7.5 across most of the star forming ring of 
NGC 7552. Thus, we conclude that velocity projection effects as¬ 
sociated with a shell-like outflow are not sufficient to explain the 
broad line profiles; an additional line broadening mechanism is nec¬ 
essary. 


4.1.2 Turbulent mixing layers 

One possible explanation is that the broad lines arise as the result 
of a turbulent gas mixing layer on the surface of cool gas clumps 
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Figure 5. Measurements of the kinematics, H a line flux, and diagnostic line ratios of NGC 7552. Vertical error bars in all figures represent Icr uncertainties, 
while horizontal error bars represent the radial bin size. The vertical dashed lines at ±0.5 kpc in all figures denote the approximate maximum radial extent 
of the star-forming ring. Measurements for all figures are taken at PA = 125°. (a): Line center velocity offsets for the narrow (blue, hats) and broad (orange) 
velocity components. The black dotted lines represent both the velocity and spatial centers of the galaxy, (b): Line widths of the narrow (blue, hats) and broad 
(orange) velocity components, (c): Star formation surface density measured in the narrow (blue, hats) and broad (orange) H a components, using the 

calibration from |Hao et al.|{201H without dust correction. The horizontal dotted line represents the threshold for driving outflows from |Heckman|{2002^ . 
(d): The ratio of broad-line H a flux to narrow-line H a flux, (e): [S II]/H a line ratios of the nan'ow (blue, hats) and broad (orange) velocity components. 
The horizontal dotted line is the threshold for shocks jHong et al.|2013t . (f): [N II]/H a line ratios of the narrow (blue, hats) and broad (orange) velocity 
components. 


l|Begelman & Fabian 1990 

Slavin et al.|1993||Westmoquette et al.| 

|2007 Tenorio-Tagle et al. 

[2010^. This theory suggests that a hot 


wind flowing out of the disk and rushing past cool gas clumps in¬ 
duces Kelvin-Helmholtz instabilities on the surface of the clumps; 
a turbulent mixing layer develops in the boundary region between 
the hot and cool gas. The broad lines would in that case be a result 
of emission from this highly turbulent boundary layer. This mech¬ 
anism has been invoked to explain the broad H a emissio n lines 
(200-300 km s”'^) seen in the central regions of M82 (e.g., West- 


|moquette et al.|2007) ). One important point is that with the addition 
of turbulent line broadening, the maximum blueshifted velocity ex¬ 
tent of the profile is no longer representative of the bulk flow veloc¬ 
ity. 


4.1.3 Shocks 

The large line widths could also be potentially explained as emis¬ 
sion from a shocked gas outflow. In this scenario, bubbles of hot 
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gas, expanding from mechanical energy injected by the starburst, 
sweep up cool, dense gas in the galaxy, propagating shock fronts 
into the ISM. The shocked shells become Rayleigh-Taylor unstable 
and fragment as they approach the scale height of the disk ( |Veilleux| 
|et al.|2005| and references therein). The observed broad line width 
in this scenario reflects the velocity of the shock ( |Ho etal.|2014) . 

Recent studies of U/LIRGs have shown widespread evidence 
of shocked outflows driven by supernova feedback from star¬ 
forming regions. These works categorize shocked emission lines by 
comparing nebular emission line ratios of [O III]/H /3, [O I]/H a, 
[N II]/H a, and [S II]/H a to shock models. [Rich et al.] ( |2011| | study 
two local LlRGs and measure broad, shocked emission lines with 
FWHM > 140 km s“^ and FWHM > 210 km s“^, respectively, 
their two targets. [Soto et ^ \2Q\2) sim ilarly measure F WHM 


Ho et al. 


2014| study 


> 235 km s“ in a sample of 39 ULIRGs. 
an isolated star-forming disk galaxy (logM, = 10.8 Mq; SFR 
~ 10 Mq yr“^) and find spatially-extended broad velocity com¬ 
ponents with line ratios that are well-explained by shock models. 
They observe a correlation between broad-component line width 
and the line ratios [N II]/H q, [S II]/H a, and [O I]/H a, and find 
that these correlations are well-explained by shock models where 
the only free parameter is the shock velocity. 

Our observations do not include the [O III], [O I], and H f5 
emission lines, but [Hong et al.j ( |2013| l calibrate shock diagnos¬ 
tics based on single line ratios. Their results reveal a dividing 
line between shock ionization and H II region photoionization at 
log([S II]/H a) > —0.5, with ratios higher than this value be¬ 
ing indicative of shock ionization. Representative examples of the 
[S II]/H a and [N II]/H a line ratios of NGC 7552 are seen in Fig- 
uresj^e) and (f), respectively. In general, the line ratios found in the 
broad emission lines are ~ 0.15 dex higher (more shock-like) than 
those seen in the narrow emission lines, but they fall slightly short 
of the threshold for shock ionization. On average, the [N Il]/Ha ra¬ 
tio increases with line with, similar to the trend found by jHo et^ 
( [2014^ . However, NGC 7552 has a much smaller dynamic range in 
line width (a = 20 — 150 km s“^) than the galaxy studied by Ho 
(ct = 10 — 350 km s“^), making detailed comparisons with models 
more difficult. 

A comparison of Figures [^c), (e), and (f), reveals a strong 
anti-correlation between H a luminosity and the [N II]/H a and 
[S II]/H a nebular line ratios. As is especially apparent in the 
[S II]/H a line diagnostic, the line ratios in the broad and narrow 
components are approximately equal in regions of very strong H a 
emission (e.g., near —0.3 kpc in Fig.j^e)). These regions are also 
where we find the most depressed [S II]/H a ratios, indicative of 
photoionization from H II regions being the dominant ionization 
source of the gas. Curiously, these regions also exhibit the strongest 
broad-line emission, with the ratio of broad to narrow H a emission 
fb/fn > 1 (Fig.j^d)). Outside of these strongest H a regions, the 
narrow-line [S II]/H a ratio continues to indicate H II region pho¬ 
toionization while the broad emission line ratio rises quickly and 
approaches the shock ionization threshold of —0.5, demonstrating 
the increased contribution of shocks in these regions away from the 
strongest photoionization. 

This anti-correlation between nebular line ratios and H a 
emission strength has been seen observationally in other LIRGs. 
[Monreal-Ibero et ^ ( |2010| > find such anti-correlations in integral- 
field spectra of 32 low-redshift LIRGs. Their analysis suggests that 
photoionization from young stars in H II regions is sufficient (or 
nearly-sufficient) in explaining this trend in only isolated, non¬ 
interacting systems. Interacting and merging systems in their study 
increasingly require the invocation of shocks, however. These re¬ 



Radius (kpc) 


Figure 6. Various measures of outflow velocity at PA = 125°. Shown are 
Au (black), the difference between the broad and narrow velocity centers, 
and the two different measurements of Aumax, Aumax,l (blue, hats) and 
Anmax ,2 (orange). The vertical dashed lines denote the approximate radial 
extent of the star-forming ring. 


suits are relevant to NGC 7552, a LIRG in the interacting Grus 
Quartet, and serve as further evidence that shocks could be a con¬ 
tributing factor to the nebular line emission we observe. 

The anti-correlation between [S II]/H a and H a luminosity in 
the context of NGC 7552 is also expected since the ages of its dom¬ 
inant nuclear star clusters are 5-6 Myr ( [Brandi et al.|2012t . The 
ionizing luminosity of an instantaneous starburst dominates the me¬ 
chanical luminosity until 8 Myr, after which the ionizing luminos¬ 
ity drops sharply while the mechanical luminosity increases until a 
few tens of Myr (Starburst99, |Leitherer et al.|1999[|Sharp & B l and^ 
|Hawthom|2010^ . The age of the star clusters in NGC 7552 places 
them within this region of dominating ionizing luminosity. Pho¬ 
toionization therefore is likely the primary driver of the nebular line 
ratios we observe, while shocks in the outflow or turbulent mixing 
layers could still be responsible for driving the broad line widths. 
Higher spatial-resolution spectra could help alleviate this ambigu¬ 
ity. [Hong et ^ j2013) note that photoionized gas is the dominant 
gas component in starburst galaxies, and caution that shock-ionized 
gas is difficult to disentangle from photoionized gas emission at the 
spatial resolution possible with seeing-limited ground-based obser¬ 
vations. 


4.2 Outflow velocities 

The net blueshift of the broad H a emission lines strongly suggests 
that the lines arise in an outflow driven by star formation feedback. 
However, it is not straightforward to infer the bulk velocity of the 
outflow from the line profile. In a simple shell model of the outflow, 
the maximum blueshifted velocity of the line profile represents the 
true de-projected velocity of the outflow. However, we cannot re¬ 
produce both the observed line widths and velocity offsets with 
such a simple model; additional line broadening from shocks or 
turbulent mixing layers appears to be necessary. With the addition 
of these line broadening mechanisms, the maximum blueshifted ve¬ 
locity is no longer indicative of the bulk outflow velocity. However, 
for consistency with the literature we proceed with the assumption 
that ttmaa: = We revisit the issue of the true velocity of the 
outflow in 

We can estimate the “maximum” outflow velocity of the gas 
based on the Gaussian line profile of the emission. Two mea¬ 
sures of Aumax are summarized by [Westmoquette et al.| ([2bl2|(: 
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Figure 7. Electron densities measured in the narrow (blue, hats) and broad 
(orange) components as a function of radius along the slit. Measurements 
from all position angles have been averaged together. The vertical dashed 
lines denote the approximate radial extent of the star-forming ring. 


A^max.l = -Abroad - A tlbroa,d(FWHM)/2 (e.g., |Rupke et al. 
|2005[ |Veilleux et al.||2005| > and Awmax ,2 = Abroad - StJbroad 
(e.g., IGenzel et al.||2011^ , where Hbroad is the velocity offset of 
the broad component line center with respect to the narrow com¬ 
ponent, Aiibroad(FWHM) IS the FWHM of the broad component, 
and (Tbroad is the Gaussian sigma value of the broad component, 
where FWHM = 2\/2 In 2 cr. We refrain from taking the absolute 
value of Aumax.i and Auniax. 2 , as in |Westmoquette et al.| ( |2012) l, 
to retain clarity about the fact that these velocities are blueshifts. 
These measures of Awmax are motivated by the fact that the emis¬ 
sion is well-characterized by a Gaussian profile. Eighty-eight per¬ 
cent of the emitting gas is therefore found at velocities red-ward 
of Aumax.i (integrating the profile from Aumax.i to oo). Sim¬ 
ilarly, 98 percent of the emitting gas is found at velocities red- 
ward of Aumax, 2 . In a spherical shell model, Aumax represents 
the deprojected velocity of the outflow. A comparison of Aumax,!; 
A?;max, 2 ; and Av, the velocity difference between the broad line 
center and the narrow line center; from PA = 125° is shown in Fig¬ 
ure]^ The velocity offset Av can be assumed to be a lower limit 
to the bulk outflow velocity of the gas. We find average values of 
(Av) = —30 ± 35 km s“^, (Awmax.i) = —180 ± 80 km s“^, and 
{Aumax,2) = —290± 100 km s“^, where the Stated 1(T Uncertainty 
is the standard deviation of the distribution of values. 


4.3 Electron density 

Our spectral window includes the [S II] 6716/6731 line pair, allow¬ 
ing us to calculate the local electron density of the emitting gas. 
This quantity is a critical component of the calculation of the mass 
outflow rate, as will be shown in |4.4| We estimate the electron 
density Ue using the [S II] 6716/6731 line flux ratio ( [OsterbroT^ 
|1989| l, calculated using the IRAF temden routine. The S/N of our 
[S II] lines is lower than that of H a and [N II], which makes an¬ 
alyzing the [S II] ratio difficult for individual spectral apertures. 
But we can still extract useful information by computing an angu¬ 
lar average of the [S II] measurements from all of our slits, similar 
to stacking the observations. By calculating this average we are 
able to compute separate electron density values for the broad and 
narrow components, as shown in Figure]^ The error bars are cal¬ 
culated by passing the allowed high and low values of the [S II] 
ratio through the temden routine. The densities measured for the 
narrow and broad components show very similar values and be¬ 


havior for |r| > 0.12 kpc. The densities of both components are 

,-3 


then drop to rie ~ 100 cm“® by |r| ~ 0.6 kpc, approximately the 
low-density limit for this diagnostic. 

These densities measured across the star-forming ring are 
higher tha n those found in typical Galactic H II regions (^ 
100 cm“^- 


Osterbrock 


1989 I. However, many starburst galaxies 


show comparable electron densities within their starburst regions. 


Annus et al. ( 1989^ measure Tie ~ 350 cm ® in 33 infrared color- 


selected LIRGs, and |Arribas et al.|j2014^ and [Soto et al.|j2012| l 
measure rie ~ 500 cm~^ in two different U/LIRG samples. In 
a study of an isolated, local disk galaxy, eTaTIpOUl mea- 
sure electron densities of a few x 100 cm“^ in a broad velocity 
component (median Tie = 300 cm“^ in the broadest component), 
with substantially lower densities in the narrow disk component 
(Tie ~ 20 cm“®). They explain the increased density in the broad 
lines by invoking compression across shock fronts, lending further 
credence to the interpretation that broad emission lines originate in 
shock environments. We do not see much difference in the densi¬ 
ties of our broad and narrow lines. This may be because the ISM in 
the high-pressure starburst region of NGC 7552 already has a high 
gas density and is not as susceptible to further shock compression. 
A stack of 14 targets in the [Newman et air| ( |2012b) sample have 
a [S II] broad-line flux ratio consistent with the low-density limit, 
and they adopt 50 cm“® for their analysis. This study is at a ~ 2 
and therefore may not be representative of the conditions found in 
starburst/outflow environments in the local universe. 

Additionally, we measure somewhat unique behavior in the 
apertures nearest the nucleus of the galaxy. Although the narrow- 
component densities drop steadily toward R — 0, which is un¬ 
surprising given the ring-like geometry of the dense star-forming 
region, the densities we measure in the broad component tend to 
peak toward the center. The broad-line measurements within ~ 65 
- 125 pc of the center show a break with the narrow-line density 
profile. The broad-component density profile shows a strong peak 
at the nucleus, although the errors on our measurements are also 
consistent with a flat profile through the nuclear region. 


4.4 Mass outflow rate 

The implications of outflows for galaxy evolution ultimately de¬ 
pend on the mass outflow rate; how much gas is being driven 
out of the galaxy due to feedback? To calculate the ionized mass 
contained within the outflow, we adopt the photo-ionized outflow 
model from jGenzel et al.]j2011[ > and jNewman et aL] ( |2Q12a| >. The 
total mass contained within the ionized outflow is given by 
1.36TTlp 

JWout = - LHa,b 

THaTle 


Mout = 6.4 X 10® M0 


10"^^ erg s 


500 cm 

Tie 



( 1 ) 


where nip is the proton mass, -/Ha is the volume emissivity of 
H a, Tie is the electron density of the outflowing gas, and 
is the H a luminosity in the broad component. Assuming a gas 
temperature of T = 10^ K, the volume emissivity of H a is 
7Ha = 3.56 X 10“^® erg cm® s“^ for Case B recombination 
|brock|198^ . With our eight different slit positions providing good 
spatial coverage in the inner ~ 1 kpc of the galaxy, we are able 
to calculate the outflow mass in a spatially-resolved fashion. We 
calculate the broad-line H a luminosity using our adopted distance 
to NGC 7552 of D = 22.5 ± 1.6 Mpc. We correct the luminosity 


Oster- 
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for Galactic foreground extinction jSchlafly & Finkbeiner||20iT| ), 
but do not make any corrections for external dust extinction, which 
we will subsequently address. Our eight slit positions provide in¬ 
complete spatial coverage, so we correct our measured Lua.b by 
the fractional spatial coverage of the observations. We also multi¬ 
ply our measured I/Ha ,6 by 2 to correct for the lack of emission 
detected from the far side of the galaxy. For the electron density Us 
we adopt the electron density corresponding to the radius of each 
spectral aperture. Due to the low S/N of the [S II] lines, this limits 
our mass outflow calculation to r < 760 pc, still large enough to 
encompass the entire projected area of the star-forming ring. Lim¬ 
iting ourselves to this radius results in 193 independent spectral 
apertures being used in our mass outflow calculation. We calculate 
separate upper and lower uncertainties for our final result by in¬ 
cluding either the upper or lower uncertainties on Ue in our error 
propagation. With these adopted values we calculate a total ionized 
gas mass of Mout = (2.5lJ |) x 10® Mq in the broad velocity 
component. 

The mass outflow rate is then the total gas mass divided by the 
outflow timescale, which we approximate as tout = 7?out/|nout|, 
the radius of the outflow divided by its velocity. The ionized mass 
outflow rate is then 

Mout = Mout/t out 


Mout = 0.2 Mq yr 


— if Mout 


VIO^Mq ) V200kms 


'i^OUt 


1 kpc 

-^out 


( 2 ) 


The radius of the outflow is a fairly uncertain measurement due 
to the galaxy’s face-on orientation. Other studies typically use the 
radial extent of broad emission as the value of Rout (e.g., |Soto| 
|et al.|201^|Rupke & Veilleux|2013| ). In the same fashion, we detect 
broad H a emission at S/N > 3 out to approximately 1 kpc in 
radius, and therefore adopt this value for flout- 

The outflow velocity is also a difficult parameter to estimate 
accurately. Numerous studies have used estimates such as Aumax.i 
and At;max ,2 as the outflow velocity, as we have done in order to be 
consistent with these studies. These measures depend on the width 
of the H a emission line, which may not be indicative of the wind 
kinematics as a whole. If the blue-shifted line centers are more in¬ 
dicative of the bulk gas velocity {Av ~ 50 km s“^), the calculated 
mass outflow rate will be smaller by factors of a few. We investi¬ 
gate the kinematics further in In order to remain consistent with 
previous studies of outflows using H a emission, we perform two 
separate calculations of Mout using two different values for Vout'- 
one using Vout = Aumax,! and one using Wout = Anmax, 2 . We fi¬ 
nally calculate mass outflow rates of Mout = 0 . 44^0 04 Mq yr“^ 
assuming Aumax,! and Mout = O.TOIq'q^ Mq yr“^ assuming 
At;max, 2 . For the sake of simplicity, we will adopt only this latter 
value for our further analysis. 

In order to represent this rate by its mass loading factor jj, 
we assume a star formation r ate in the range of M* = 10 - 
15 Mq yr“^, the estimate from Pan et al. (2013 ' of the SFR of the 


starburst ring. Assuming the outflow is driven primarily by feed¬ 
back from the central starburst, it is reasonable to adopt the SFR of 
that region for our calculations, rather than an estimate of the SFR 
from the whole disk. Our estimate of the mass outflow rate corre¬ 
sponds to mass loading factors in the range of 77 = 0.05 - 0.07 with 
no extinction correction. 

The |Arribas et al.| p014| l study finds average mass loading 
factors of Tf = 0.3 and 77 = 0.5, respectively, for the LIRGs and 


ULIRGs in their sample, with very few systems having 77 > 1. 
These averages are factors of a few higher than our measurements 
in NGC 7552 but have been corrected for extinction, unlike our 
data. Our spectra do not include H /3, so we are unable to do a di¬ 
rect calculation of the broad component extinction. Without a direct 
measurement, correcting our data for extinction becomes a much 
more imprecise task. 

Direct measurements of the extinction toward embedded star 
clusters in the star-forming ring of NGC 7552 are presented by 
[Brandi et al.|P01^ . The most embedded clusters have measured 
extinctions as high as Ay = 9.0, with an average value of Ay = 
5.3 over the star-forming ring. However, this extinction estimate 
represents the attenuation by the dense gas within the star-forming 
ring, along a line-of-sight extending all the way to the deeply- 
embedded star clusters. Emission from supernova-driven outflows 
on the near side of the disk would not suffer from such strong atten¬ 
uation, and these extinction values are therefore likely to be over¬ 
estimates. 

We also have an estimate of the extinction from our fitting of 
the UV stellar continuum (j ]3.2^ . Assuming a jCalzetti et al.|(2000l > 
extinction law, this extinction of E{B — V) = 0.64 attenuates our 
H a flux by a factor of 7.1. That the star cluster along the COS line- 
of-sight is UV-bright implies that it has broken out of the dense gas 
in the galaxy’s mid-plane and star-forming ring, and the attenuation 
of the starlight is largely the result of extra-planar gas (e.g., gas 
entrained in a galactic wind). Adopting this extinction results in an 
adjusted mass outflow rate of Mout = ^ Mq yr“^, or a mass 

loading factor of 77 = 0.3 - 0.5, similar to the jArribas et al.|p014| l 
study. We adopt this result as our best estimate of the mass outflow 
rate. 


5 MEASURING THE OUTFLOW USING UV 
ABSORPTION LINES 

One concern with our analysis is that the gas tracers we have used 
are only measured in emission, a process that depends on the square 
of the local gas density. This has the effect of making our study 
more sensitive to denser regions of the outflow, and not necessarily 
representative of the outflow on a large scale. For this reason it is 
informative to compare our results to UV absorption-line measure¬ 
ments of NGC 7552. A study of the mass outflow rate using this 
data was also presented in jLeitherer et al.H2013^ . A major differ¬ 
ence between that study and our work is that they estimated a total 
gas column from the measured dust attenuation. In this study, we 
pursue an alternative method, using column densities derived from 
the absorption profiles of Si lines in order to calculate a total mass 
outflow rate. 


5.1 UV absorption line profiles: gas at 1000 km s ^ 

In Figure we show the HST/COS spectrum of NGC 7552. The 
spectrum shows many prominent blueshifted interstellar absorp¬ 
tion lines. To study these lines, we first normalize the spectrum by 
a theoretical model of the stellar continuum as described in §3.2| 
The spectrum and fitted continuum model are shown in the top 
panel of Figure This continuum model accounts for broad ab¬ 
sorption from stellar winds that could contaminate our target inter¬ 
stellar lines. The lower left panel of Fig. shows the line profiles 
of Ly a, C II 1335, and Si II 1260. These are some of the strongest 
and most saturated transitions in this wavelength range and they 
are therefore useful for exploring the full velocity extent of the gas 


© 2015 RAS, MNRAS OOO.pTfT^ 


























Supernova-Driven Outflows in NGC 7552 11 




-1000 -500 0 500 -1000 -500 0 500 

Velocity (km s“') Velocity (km s“^) 


Figure 8. HST/COS spectrum of NGC 7552. The top panel shows the data (black) and a theoretical stellar population model (blue) from STARBURST99. 
Prominent ISM absorption lines are labeled. Emission lines other than the marked Ly a line are geocoronal lines. The lower left panel highlights the continuum- 
normalized spectra of Ly a, C II 1335, and Si II 1260. These features are some of the strongest and most-saturated transitions in the spectrum, and they show 
the full extent of the absorption, extending to > 1000 km s“^. The absorption profile is very deep (almost reaching zero intensity) from 0 to —600 km s“^, 
suggesting that the gas in the outflow almost completely covers the source. In the lower right panel, transitions of Si II, Si III, and Si IV are shown. The 
similarity of the line profiles suggests that different ionization phases of the gas have similar spatial distributions and kinematics. 


associated with the wind. Absorption is clearly seen in all ions out 
to 900 km s“^ and in Ly a out to 1100 km s“^. These velocities 
greatly exceed the maximum velocities associated with the H a 
emission (Atimaa; = 160 —250 km s“^), a point that we will return 
to in The absorption is very deep from ti = 0 to —600 km s“^, 
almost reaching zero intensity. This implies that the gas in the wind 
at these velocities almost completely covers the UV-luminous stars. 

The Ly a line shows a relatively typical profile for a star form¬ 
ing galaxy: redshifted emission and blueshifted absorption. Such a 
profile can arise when Ly a photons resonantly scatter in the wind. 
The redshifted emission is thought to be produced by photons scat¬ 
tering off the opposite side of the wind (e.g., the portion of the 
outflow moving away from the observer). It is somewhat surprising 
that we detect Ly a. emission in our COS observations of NGC 7552 
since we are observing a UV-bright cluster complex against a back¬ 
drop of very dusty, heavily-embedded clusters. [Brandi et al.H2012) 
compute an average dust attenuation for NGC 7552’s star forming 
ring of Av = 5.3, which corresponds to attenuation by more than 
a factor of a million at 1216 A. Since the ring is effectively opaque 


in the UV, it is unlikely that the Ly a. emission comes from the far 
side of a galaxy-scale outflow. It is more plausible that the emis¬ 
sion arises from the back side of a cluster-scale outflow (e.g., from 
gas on the near side of the disk that is being driven away from the 
observer and towards the disk mid-plane.) |Leitherer et al.| [2013| > 
model the Ly a profile of NGC 7552 using the radiative transfer 
models of |Schaerer et al.H2011[ . They find a good fit to the emis¬ 
sion profile using a spherical shell model with a constant expansion 
velocity of 50 km s“^, which is comparable to the velocity of the 
broad H a component. However, the model substantially under¬ 
predicts the high velocity Ly a absorption. One plausible explana¬ 
tion is that the absorption traces both a low velocity cluster-scale 
outflow and a higher velocity galactic-scale outflow. We consider 
this idea further in section 


5.2 Measuring column densities 

Our goal is to use the UV lines to estimate the column density of 
the outflowing gas and, ultimately, to compute the mass outflow 


© 2015 RAS, MNRAS 000,[T|fT^ 






















































12 Corey M. Wood et al. 



1000 

-500 0 

-1000 

-500 


Velocity (km s“') 


Velocity (km s“^) 


Figure 9. Line and column density profiles derived from various Si If and S IV transitions. The transitions and their oscillator strengths are listed in the legend. 
Portions of the line profile that are contaminated by other lines are shown as dotted lines. The disagreement between the derived column densities of Si II may 
be due to a non-unity covering factor or the presence of unresolved saturated components. The low /-value transition of Si II 1304 is likely lo yield Ihe most 
reliable column density, however the blue portion of the profile is contaminated by 0 I 1302. The dashed black line shows the column density derived from 
Si II 1260 arbitrarily scaled up by a factor of 8. Since this provides a good match to the column derived from Si II 1304, we use this to estimate the total Si II 
column. The Si IV lines are not strongly saturated and yield column densities that are in reasonable agreement. 


rate. The first step is to derive accurate ionic column densities for 
all of the important Si II phases and combine them to arrive at a 
total Si column. By applying a correction for depletion onto dust 
and assuming a solar Si/H ratio, we can estimate the total hydrogen 
column. The right-hand panel of Fig. shows transitions of Si II, 
Si III, and Si IV, which have ionization potentials of 8.15, 16.35, 
and 33.49 eV, respectively. The overall shapes of the line profiles 
are very similar, suggesting that all of the lines are associated with 
the warm neutral phase of the ISM. Unfortunately, the strongly sat¬ 
urated nature of the Si II 1260 and Si III 1206 lines makes it difficult 
to derive reliable ionic column densities. 

In Figure]^ we show column densities derived using the ap¬ 
parent optical depth method jSavage & Sembach||199l| > for the 
Si II 1190, Si II 1260, and Si II 1304 lines, which have /-values that 
differ by as much as a factor of 14. In the ideal case, all three tran¬ 
sitions would yield the same column density. The fact that the lines 
with lower /-values have higher derived column densities may in¬ 
dicate that unresolved saturated components are present jSava g e &| 
|Sembach|1991^ or that the covering factor of the gas is less than 
unity. The low /-value transition of Si II 1304 is likely to yield the 
most reliable column density, however the blue portion of the pro¬ 
file is contaminated by O I 1302 geocoronal emission. We find that 


scaling up the column density derived from the strong unblended 
Si II 1260 line by a factor of 8 provides a reasonably good match to 
the column derived from Si II 1304, and we use this rescaled line 
profile to estimate the total Si II column. The sole Si III transition 
A1206 has a high oscillator strength and is likely strongly saturated. 
The Si IV 1393, 1402 lines do not appear to be strongly saturated 
and yield column densities that are in reasonable agreement with 
one another. 

The ratio of Si II to Si IV is approximately 5:1. Assuming that 
the gas arises from a single phase, we can use the Si Il/Si IV ra¬ 
tio to estimate the contribution of Si III to the total Si column. We 
compare our measured ratio to the photoionization model predic¬ 
tions for circum-galactic gas presented in |Werk et al.|pOT4) . We 
use the model shown in their Fig. 41 which is for log N(H I) = 
17. In the model, the relative ratios of different ions depend upon 
the ionization parameter. We constrain the ionization parameter us¬ 
ing our measured ratio of Si Il/Si IV and infer N(Si III) = 2x 
N(Si II) from this model. This photoionization model was designed 
for circum-galactic gas illuminated by the diffuse UV background, 
and therefore may not be an accurate representation of the gas in 
the inner parts of the galactic wind, which is likely a complex, 
multi-phase structure. Additionally, this model assumes a metal- 
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licity of log{Z/ZQ) = —0.4, likely lower than the metallicity of 
NGC 7552. 

We compute the total Si II column as a function of velocity 
by summing the Si II, Si III, and Si IV profiles and find N(Si) = 
9.86 X 10^® cm~^. We turn the Si column into an H I column 
by applying a correction for dust depletion and assuming a solar 
Si/H ratio, using values for these parameters of D(Si) = —0.29 
and log(Si/H )0 = —4.45 from|Draine|(|201 l|l. The total hydrogen 
column density in the outflow, estimated from the Si II, Si III, and 
Si IV absorption lines between velocities of 0 and —1000 km s“^, 
is N{H) = 5.28 X 10^° cm'^ 

One disadvantage of the apparent optical depth method is the 
implicit assumption that the gas completely covers the source. In 
principle, one can obtain information about both the covering fac¬ 
tor and the column density by simultaneously analyzing multiple 
transitions of a given ion that span a range of /-values. We ex¬ 
plore this approach by simultaneously fitting the Si II 1190,1193; 
Si II 1260; and Si II 1304 lines with Voigt profiles. The observed 
line profiles are complex and clearly contain multiple blended com¬ 
ponents. Given our modest spectral resolution (~ 100 km s“^), we 
model each absorption feature with a set of seven individual veloc¬ 
ity components. Each velocity component is fit with a Voigt profile 
with the fe-parameter set to half the instrumental resolution. The 
covering factor and Si II column density are free parameters for 
each component. 

The similarity of the Si II, Si III, and Si IV line profiles (see 
Fig. motivated us to include the Si III 1206, and Si IV 1393, 
1402 lines in our fits. We allow the ionic ratios (e.g.. Si IV/Si II) 
to vary, but we assume that this ratio is the same for all velocity 
components. With these assumptions, we achieve excellent fits to 
all of the individual absorption line profiles. Once again, the Si III 
column is poorly constrained; we therefore assume N(Si III) = 2x 
N(Si II) as we did previously. Summing over all blueshifted com¬ 


ponents we findA''(Si) = (1.71 ±0.45) x 10^ 


' which yields 


N{H) = (9.17±2.41) X10^° cm“^ using the same abundance and 
dust depletion corrections as above. This column is roughly a fac¬ 
tor of two larger than we derived using the apparent optical depth 
method, but is consistent with that measurement within the errors. 
Most of this increase comes from a much larger derived column at 
low velocities (< 200 km s“^). 

[Leitherer et al.| ( |2013^ estimate the hydrogen column density 
in a different manner, using the galactic relationship between N(H) 
and reddening from |Diplas & Savage| ( |199^ 

N(H) = 4.9 X 10^^ - V) cm"^ (3) 

Using the E (B — V) = 0.64 derived from fitting the stellar contin¬ 
uum (1 3.2 1 , we calculate N(H) = 3.1 x 10^^ cm“^, about a factor 
of 5.9 and 3.4 higher than the value derived from the Si absorption 
lines using the two methods described above. It is reasonable to ex¬ 
pect that this value represents an upper limit on the column density, 
as the extinction traces the column to the star cluster and thus likely 
includes a contribution from both gas in the disk and the wind. To 
be maximally conservative, we adopt the lowest of our three H I 
column density estimates, N{H) = 5.28 x 10^° cm“^, derived 
using the apparent optical depth method. 


5.3 The wind mass and mass outflow rate 

Following |Rupke et al.| p005| >, we model the outflow as a thick 
shell with inner and outer radii ri and r 2 . The total mass of the 
outflow is computed as 

t^wind — 4:TvCnCffirnpN{H)rir2. (4) 



-1000 -500 0 

Velocity (km s”') 



Figure 10. Outflow model fit to NGC 7552’s UV absorption lines. The 
model, from |Steidel et al!]j2QlQ) , requires both the cloud covering factor 
and acceleration to have a power law dependence on the distance from the 
starburst. In the top panel, the gray and black lines show the Si II 1260 and 
C II 1335 absorption lines and the red line is the best fit model. The bot¬ 
tom panel shows the resulting dependence of the velocity (solid) and cloud 
covering factor (dashed) on radius. Note that most of the gas acceleration 
happens over a very small range of radii (1-3 kpc). The minimum and 
maximum radius (1 and 100 kpc) are not free parameters of the fit. 


Here Cn represents the angular covering fraction of the wind, C/ 
the clump covering factor, p the mean molecular weight, and nip 
the proton mass. We adopt Cn ~ 0.67, assuming that the wind 
has a bi-conical geometry with a half-opening angle of 60° ( |Chen| 
|et al.|201^ . (Note that |Rupke et al.|p005^ assume a slightly smaller 
value of Cn ~ 0.4 for LIRGs). The clump covering factor C/ 
accounts for the fact that the wind contains discreet clouds that do 
not completely cover the UV emitting stars. Based on our multi- 
component fits to the line profiles we find that the covering factor 
ranges from 0.6 - 1. We adopt the average value, 67/ = 0.9, in our 
calculation. We assume — 1.4 to account for both hydrogen and 
helium. The outflow radii are the least well-constrained parameters 
of the equation. For the inner radius we adopt ri = 1 kpc, the 
radius we assumed for our H a emission in §4.4| Given that the 
velocity extent of the UV absorption is much larger than that of the 
H a emission, it is reasonable to assume that the outer radius of the 
shell is substantially larger than that traced by H a.|Chisholm et al.| 
( |2014^ adopt a wind radius of 4.4 kpc based on an observed scaling 
between star formation and outflow extent ( [Grimes et al.||2005l ). 
We adopt this value as our outer radius r 2 .|Rupke et al. (|2005|l and 
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other authors, including |Leitherer et al.| ( [2013| l, frequently assume 
r = 5 kpc, although these particular studies use thin shell models 
(n « r 2 ). Rewriting the above formula in more convenient units 
and plugging in values we find 


= 2.0 X 10* MfT) f — I ( — I X 
^ V 0.67/V 0-90y 

N{H) \ / n \ / r2 

5.28 X 102°cm~2/ ^Ikpc/ \ 4.4 kpc 

To compute the average mass outflow rate, we divide the wind mass 
by the average timescale of the wind, t = r 2 /v, where v is the aver¬ 
age outflow velocity and r 2 is the outer shell radius. Here we have 
adopted the column-density-weighted average of n = 410 km s“^, 
which yields a timescale of 10 Myr. Plugging in we find 



M_. = 19M0yr- 

N{H) 


0.90 


5.28 X 1020 cm-2 


ri 

1 kpc 


V 


410 km s 


( 6 ) 


implying a mass loading factor in the range 77 = 1.3-1.9. This con¬ 
trasts with the |Leitherer et al!]j2013^ study, which estimates mass 
outflow rates on the order of ~ 100 Mq yr~^ for their sample. 
Their result could be an over-estimate, given that the dust attenua¬ 
tion is measured in the starlight, while the high-velocity gas likely 
exists at larger disk heights. There is also significant uncertainty 
in our measurements given the large ionization corrections for the 
various Si species. 

In practice, it is likely that the outflow spans a range of radii. 
ISteidel et al.| ( |2010l l develop a model of the kinematics and geom¬ 
etry of circumgalactic gas based on extensive observations of UV 
absorption lines in projected galaxy pairs ai z = 2 — 3. In their 
model, cool dense clouds are embedded in a hot, diffuse, volume¬ 
filling medium. These clouds are accelerated over time by radiation 
pressure and ram pressure from the hot wind. As they move out¬ 
ward, the external pressure supplied by the hot wind declines and 
the clouds expand. But they do not expand by enough to counteract 
their geometrical dilution, and thus the cloud covering factor drops 
with radius as Cf oc r \ with 7 = 0.2 —0.6. For the Lyman Break 
galaxies modeled by Steidel, typical maximum values of the veloc¬ 
ity and radius are Vmax ~ 750 km s“^ and r^ax ~ 80 kpc. This 
outer wind radius is significantly larger than the radius we have as¬ 
sumed (r = 4.4 kpc). We have experimented with fitting Steidel’s 
analytic model to the absorption profile of NGC 7552 in order to 
determine the relationship between velocity and radius. However, 
we were unable to produce a satisfactory fit due to the fact that the 
cloud covering factor in NGC 7552 is roughly constant over a h uge 
range of velocity (n = 0 to —600 km s“^). As shown in Figure 10 


the best fitting model requires a very rapid acceleration of the gas 
over a very small range of radius (1-3 kpc), which may be difficult 
to achieve. Unfortunately, the model has many free parameters and 
it is not possible to uniquely constrain the inner and outer radius of 
the outflow (e.g., equally good fits can be obtained by setting the 
starting radius to 0.1 kpc or the outer radius to 10 kpc). What is 
clear is that it is a challenge to devise a physical model to explain 
how the cloud covering factor remains roughly constant over a very 
large range of velocity. 


6 DISCUSSION 

We have measured the outflow mass and mass outflow rate of the 
nuclear starburst region of NGC 7552 using both H a emission and 
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Figure 11. Comparison of the UV absorption and H a emission showing 
the ionized silicon column density inferred from Si II, Si III, and Si IV UV 
absorption lines (periwinkle), the coincident H a profile (black dotted), and 
the fit to the broad H a velocity component (red). The H a data has been 
truncated at 500 km s"'^ due to contamination from the [N II] 6548 line. It 
is important to note that the Si absorption traces only gas in front of the star 
cluster along our line-of-sight, while H q is sensitive to gas both in front 
and behind the cluster, largely dependent on only extinction effects. 


UV absorption lines. For the H a line, we measure the mass outflow 
by decomposing the visible spectrum into broad and narrow veloc¬ 
ity components. We detect broad-line emission along all sightlines 
within ~ 1 kpc of the galactic nucleus. The narrow-line emission 
shows a reasonably well-behaved rotation profile. The broad emis¬ 
sion also largely traces this rotation pattern, but is blueshifted by 
Aw ~ 30 - 50 km s“^. We find that the broad component has line 
widths of 240 - 350 km s“^ FWHM, while the line widths in the 
narrow component are ~100 km s“^. Our measurements of the 
narrow component are consistent with it arising from star-forming 
H II regions within the disk. We attribute the broad component to a 
feedback-driven galactic wind, similar to recent studies (e.g., |Ne\^ 
[man et al.| 2012 b| l, and calculate an ionized mass outflow rate from 
the broad H a flux. Combined with a measurement of the local 
electron density via the [S II] line ratio diagnostic, this mea¬ 
surement yields a mass outflow rate of Mout = 5 . 02 o '5 M 0 
with a mass loading factor of ?) = 0.3 - 0.5. 

Our measurements show electron densities of Tie ^ 500 cm“* 
in the broad lines, substantially higher than what has been mea¬ 
sured or assumed in some recent studies (typically Tie ~ 10 - 


20131. We 


50 cm~*) I Newman et al.|[^12b Rupke & Veilleux 
note, h owever, that the measurement of rie by [Newman et al.| 
1 2012b I (rie = 102io° crn~*) is still statistically consistent with 


our measurements, considering the large uncertainty in their value. 
That study also measures Tie using a spatially-integrated stack of 
14 galaxies, and the resulting [S II] ratios may be sensitive to re¬ 
gions outside the main starburst of these systems. Many local stud¬ 
ies of (U)LIRG/starburst systems find elevated electron d ensities 
(Tie ^ 300 cm“*) in the central star-forming regions (e.g.,|Armus 


|et al.|198^ [Lehnert & Heckman||1996[ [Weistrop et al.|201^ Ar- 
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[ribas et al.|2014| ), as well as large declines in tie with radius. These 
high central densities and radial dependences suggest measuring Us 
in a spatially-resolved fashion may be important for constraining 
the densities of the starburst region. Given that Eqn.[2 depends in¬ 
versely on He, increasing the precision in measuring rie could have 
order-of-magnitude effects on the measurement of Mout. This may 
explain the largest source of the discrepancy between our measure¬ 
ments and those at high-redshift, although it is not clear that winds 
at high redshift are analogous to local starburst-driven winds. 

A main source of uncertainty in our calculation of Mout is ex¬ 
tinction in our observations because our spectra contain no reliable 
method for estimating the extinction. Our spectral window does not 
include H j3 and we are therefore unable to directly measure the ex¬ 
tinction via the Balmer decrement H a/H /3. Our best estimate of 
the extinction in the broad H a component is made by using the red¬ 
dening of starlight measured via our UV stellar continuum model 
fitting (jp^. 

We have also analyzed a far-UV HST/COS spectrum of 
NGC 7552 and measured properties of the outflow via absorption 
profiles. We converted absorption profiles of Si II and Si IV to a 
total hydrogen column density and recovered a total gas mass by 
assuming a simple radial shell geometry. This yields a total mass 
contained within the wind of M^ind = 2.0 x 10® Mq with a mass 
outflow rate of Mmind = 19 Mq yr“^, implying a mass loading 
factor of = 1.3 - 1.9. We use a conservative estimate of the col¬ 
umn density for this calculation; the true column — and, similarly, 
its associated mass outflow rate — could be 2 to 6 times higher. 
This mass outflow rate is a factor of ~ 2 higher than our mass out¬ 
flow rate measured with the optical emission lines, once that value 
is corrected for extinction. 

Compared directly, these two different analyses show that the 
UV and H a emission have very different kinematic signatures and 
trace different regimes of the outflow. The UV absorption profiles 
extend to velocities upward of —1000 km s“^, far higher than seen 
in the H a profile. Our most liberal calculation of the maximum 
H a outflow velocity has an average value of only {Aumax, 2 ) = 
— 290 ± 100 km s“^. As seen in Figure [IT| the majority of the gas 
column traced by the UV absorption exists at velocities larger than 
this. There exists a large amount of material at high velocities that 
our H a observations simply do not detect. The ionization state 
of the gas inferred from the Si Il/Si IV ratio suggests that > 99 
percent of the hydrogen is ionized. The lack of H a emission at 
these velocities is therefore likely due to a large drop in the local 
gas density, rather than a dearth of ionized hydrogen. In short, H a 
appears to be insensitive to a large fraction of the outflow mass of 
a similar gas phase. 

It is somewhat surprising that we see strong, redshifted Ly a 
emission along our line-of-sight (Fig.[^ lower-left). In a galactic 
wind scenario, redshifted emission would be the result of backscat- 
tering from outflowing material on the far side of the galaxy. But 
the disk, with an average extinction of Av = 5.3, attenuates light 
by more than factor of a million at 1216 A, and thus no Ly a emis¬ 
sion from the far side of the disk would be present in our observa¬ 
tions. One logical explanation of this emission is that it traces ma¬ 
terial outflowing at a star-cluster scale, rather than tracing a bulk, 
galactic-scale outflow. [Leitherer et al.jpofs) model the Ly a profile 
as a spherically-expanding shell with a constant expansion velocity. 
They find a reasonable match to the emission profile using an ex¬ 
pansion velocity of 50 km s“^. The model fails to account for sub¬ 
stantial absorption at high velocities, however (see their figure 18). 
One concern is that this excess absorption is the result of Galac¬ 
tic Ly a absorption contaminating the line profile. However, even 


the metal lines are nearly saturated beyond —600 km s“^, where 
their model accounts for only ~ 50 percent of the Ly a absorp¬ 
tion. Therefore, the shell model and Galactic Ly a contamination 
must not account for all the absorption we see at high velocities. 
It is plausible, then, that the Ly a simultaneously traces multiple 
outflow components: a low-velocity (~ 50 km s“^) cluster-scale 
outflow and a higher-velocity (up to 1000 km s“^) galactic-scale 
wind. 

Our interpretation of the H a emission is that it likely arises 
on the cluster scale, similar to the Ly a emission, rather than being 
associated with any global galactic wind. In addition to the large 
disparities in the velocity profiles of H a and various UV tracers, 
this interpretation is supported by the high electron densities seen 
in the broad emission lines. Most of the broad H a emission within 
the starburst ring comes from gas with electron densities of 300 
- 700 cm“®, factors of a few higher than typically found in H II 
regions but consistent with other starbursts. It is unlikely that any 
gas clumps entrained in a wind would maintain these high densities 
to large radii. jWerk et al.|(|2014J infer typical gas densities for cool 
{T = K)CGMofne = 10“® - lO"'^ cm"® at distances of 20- 
150 kpc from the host galaxy. It is more likely that these gas clumps 
exist near the disk and are being swept up in shocked supemovae 
outflows. Based on the common picture of galactic winds — hot 
gas mshing at high velocity past dense, warm clouds of gas — the 
optical emission lines are likely arising from warm clouds that have 
not yet been accelerated by the wind. 

Looking at the H a emission in isolation, these clouds could 
either be material localized to star-forming regions and not yet en¬ 
trained in the wind, or newly-entrained material at the base of the 
wind that has not yet been dismpted and accelerated to high ve¬ 
locities. The broadening of the emission lines is then due to either 
turbulent mixing on the surface of the clumps (Begelman & Fabi^ 


|1990[[Slavin et al.|1993[|Tenorio-Tagle et al.|2010^ or due to shock 

processes at the boundary of a local supernova bubble. The line 
ratio diagnostics (Figs. |^e),(f)) demonstrate that the broad lines 
suggest shock-like ionization over much of the starburst ring. The 
[S II]/H a line ratio, a more reliable single-line diagnostic, sug¬ 
gests that shocks could play a role in the gas ionization, with ra¬ 
tios near the shock ionization threshold (log([S II]/H a) ~ —0.5) 
in all but the strongest regions of H a flux. Shock signatures can 
be diluted by strong photoionization in ground-based observations 
( [Hong et al.|2013f , and our line ratio measurements are therefore 
likely to be underestimates. 

Although we do not have information about how high above 
the disk this broad-line gas is located, we look at the nearby edge- 
on starburst galaxy M82 to shed light on how the H a-emitting gas 
is distributed as a function of disk height in a starburst galaxy. In 
Figure [T^ we show the H a profile of M82 as a function of height 
above the mid-plane, summed in a 1 kpc-wide box centered about 
the nucleus. We use a narrow-band H a image of M82 acquired 
with the Mini-Mosaic Imager on the WIYN 3.5 m telescope ( |West-| 
|moquette|2007^ . We sum only the flux in the southern outflow, as 
the gas near the mid-plane in the northern outflow is partially ob¬ 
scured by the disk. This narrow-band imagery does not distinguish 
between the narrow or broad components, but we find that the broad 
emission in NGC 7552 typically accounts for 50 percent of the total 
H a flux. If we assume this is also the case in M82, we can attribute 
half of the flux to the narrow emission. If we further assume that 
the narrow emission arises in the disk, and thus all emission at low 
disk height is in the narrow component, then all narrow emission 
(50 percent of the total) is located within 160 pc of the mid-plane 
in M82. Gas within 1 kpc accounts for 94 percent of the total H a 
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Figure 12. Left: Log-scaled (solid) and cumulative fraction (dashed) of H a counts in M82, as a function of projected height above the disk. Counts are 
summed within a 1 kpc-wide box centered about the nucleus. Right: Continuum-subtracted H a grayscale image of M82 jWestmoquette|2007^ . Image is 
oriented with north up and east to the left. The orange box represents the 1x2.5 kpc extraction region. 


emission, or 87 percent of the remaining emission after discount¬ 
ing the first 50 percent attributed to the narrow component. This 
supports the idea that a majority of the gas in the broad velocity 
component is located near the disk in starburst outflows. This is 
also consistent with the high densities in the broad component, as 
it is unlikely that the T = 10^ K gas would remain at such high 
densities out to large disk heights, and we therefore expect to pref¬ 
erentially detect gas near the disk. 

The velocity centers of the broad component also support the 
idea that this component exists near the disk. Comparing the veloc¬ 
ity centers between the broad and narrow components in Fig.|^a) 
shows that the broad components still largely trace the rotation 
curve of the galaxy, despite their modest blueshifts. Angular mo¬ 
mentum conservation would be expected to substantially lower the 
rotational velocity of any outflowing gas as it reaches large radii. 
The fact that the broad component still shows strong signs of ro¬ 
tation therefore suggests that the emitting gas remains close to the 
disk. 

We have additionally created a toy model of an expanding 
spherical shell of gas, and compare the resulting velocity profile 
to our typical broad-component H a line profile. Coded in IDL, the 
model creates a three-dimensional cube of cells, with cells falling 
within defined inner and outer radii considered to have an arbi¬ 
trary unit of gas in them. We compute the column density for each 
cell by summing all the gas from that cell toward the observer. We 
then normalize the column through the center of the sphere to the 
measured UV reddening E{B — V) = 0.64, and calculate the 
extinction toward each cell using a |Cardelli et al.| ( [l989^ model. 
Each cell has a constant radial velocity that is projected into the 
line-of-sight, and then we create a velocity profile by construct¬ 
ing a histogram of the projected velocities. We convolve this pro¬ 
file with a Gaussian profile representing the combined effects of 
our instrumental resolution (a = 17 km s“^) and turbulent/shock 
broadening. We find that a constant-velocity shell with an expan¬ 
sion velocity of 100 km s“^ that is turbulently-Zshock-broadened 
by (T = 90 km s“^ creates a line-of-sight Gaussian velocity profile 
with (j = 102 km s“^ (~ 240 km s“^ FWHM) and a blueshifted 


centroid of Av = —48 km s“^. These values agree remarkably 
with our broad component measurements. 

This model shows that the velocity center depends strongly on 
E(B — V), and that the line width depends strongly on local broad¬ 
ening effects. We have experimented with somewhat more-complex 
models by varying the radial density and velocity structures and 
find broadly similar results. Varying the outflow velocity in the 
model leads to roughly linear changes in the velocity centroid off¬ 
set, and an outflow velocity of 200 km s“^ has a profile blueshifted 
by over 100 km s“^, much larger than what is seen in our data. 
We also experimented with adding an additional velocity compo¬ 
nent by including a second shell with a different constant outflow 
velocity. A second shell with a lower outflow velocity in our model 
creates a combined profile with a lower blueshifted centroid veloc¬ 
ity, but that is still largely consistent with our data. If the additional 
shell has a larger outflow velocity, the velocity centroid increases 
rapidly in response: a second shell with Wout = 200 km s“^ cre¬ 
ates a combined profile with a centroid of —87 km s“^, inconsistent 
with our data. Our model therefore suggests the line profile largely 
reflects the fastest outflow component, ultimately resulting in an 
over-estimate of the average outflow velocity of the gas. 

While this model may not be a perfect representation of the 
data, it demonstrates that simple scenarios involving low-velocity 
{v ~ 100 km s“^) expanding gas shells can plausibly recreate 
the observed broad-component H a kinematics with contributions 
from local broadening processes. This model may not represent 
the exact physical scenario present in the star-forming region of 
NGC 7552, but we believe it demonstrates that such low-velocity 
scenarios need to be considered before adopting outflow veloci¬ 
ties factors of a few larger than this model value. If we adopt 
t his model expansion velocity for the outflow velocity Vout in 


[4A our calculated Mout is reduced by nearly a factor of 3 to 
Mout = 2.2 Mq yr~^ (tj — 0.1 - 0.2), a factor of 8.5 smaller 
than the UV measurement. We believe this Vout represents a more 
accurate measurement of the bulk velocity of the gas, and there¬ 
fore we adopt this adjusted Mout as our best estimate of the mass 
outflow rate. 


© 2015 RAS, MNRAS OOO.pTfT^ 











Supernova-Driven Outflows in NGC 7552 17 



Figure 13. Schematic diagram of our inteipretation of the broad nebu- 
iar emission in NGC 7552. Star formation in cioseiy-spaced ciusters in the 
disk create hot, pressurized hubbies due to strong winds and supemovae. 
These hubbies can break out of the disk, iaunching dense materiai outward. 
These dense gas ciumps expand as they travei to iarger distances above the 
disk, iowering the iocai density and quickiy diminishing the amount of H a 
being emitted. The iine widths of the nebuiar iines are iikeiy the resuit of 
unresoived shocks or turbuient motions within these ciumps, rather than 
representing the projection of the buik flow veiocity of the ciumps. 

If the line widths of the broad H a velocity component are 
being driven largely by surface layer turbulence or shock pro¬ 
cesses, then the line profile would contain little useful information 
about bulk motions of the gas. Other recent observational work 
has reached similar conclusions about the mechanism behind the 
broad line widths. |Ho et ah] ( |2014^ present an integral field study 
of a single isolated disk galaxy with similar global properties to 
NGC 7552. The properties of the broad H a components in their 
study largely resemble those which we find in NGC 7552, e.g., 
large line widths, high electron densities, and shocked line diag¬ 
nostics. They reach the conclusion that the line ratios and velocity 
dispersions of the broad velocity component are well-represented 
by shock models, and that compression along shock fronts leads to 
the high electron densities as measured via the [S It] line ratio. 

We show a cartoon schematic of our physical picture in Fig- 
ure[T^ Vigorous star formation is on-going in closely-spaced clus¬ 
ters in the disk, and strong stellar winds and supernovae create 
hot, pressurized bubbles within the disk, expanding at ~ 50 — 
100 km s“^. These bubbles can break out of the top and bottom 
of the disk, launching dense material outward. This gas is likely 
shocked, which drives the nebular line ratios and compresses the 
gas, leading to higher local gas densities compared to the star form¬ 
ing regions from which they originate. These dense gas clumps are 
likely at high pressures (rie ~ 500 cm~^ and T ~ 10^ K implies 
pressures of P/k > 10®), and would then expand as they travel 
to larger distances above the disk, lowering the local density and 
quickly diminishing the amount of H a being emitted. 

Recent simulations of galactic winds have also studied the 
warm gas in supernova-driven winds. [Creasey et al.| ( |2013) l model 
a simplified ISM with star formation following the Kennicutt- 
Schmidt relation, with frequent supernova events in the disk driv¬ 
ing gas into the halo. They find that nearly all of the warm gas 


(T = 10^ K) in their simulation exists in dense clumps at low 
velocities (v < 350 km s“^). They go on to note that the interface 
between the warm gas and the hot wind in their simulation could be 
the source of lower ionization lines like [S II]. [Cooper et al.|f2009[ l 
model the evolution of dense, radiatively-cooling gas clouds within 
a hot wind. Their results show such a cloud breaks into numer¬ 
ous small, dense “cloudlets” due to Kelvin-Helmholtz instabilities 
at the boundary between the dense gas and the hot wind. These 
cloudlets typically exist at outflow velocities of n ~ 200 km s“^ 
in their simulation. However, their emission-weighted 2-velocity 
profile of the H a line along a filament shows the emission peak¬ 
ing at 15 km s“^. This may be qualitatively similar to the picture 
painted by other H a studies: the “outflow velocity” of the H a gas 
is substantially larger than the velocity of the H a profile centroid. 
Neither study shows the presence of warm gas at high velocities, 
however, as seen in our UV data. But neither study follows the gas 
out to large radii. 

Overall, our results simply suggest that caution must be ex¬ 
ercised when interpreting broad H a velocity components as large- 
scale galactic winds. The broad H a component in NGC 7552 likely 
arises from dense material either being shocked or undergoing tur¬ 
bulent mixing. It is still plausible — even likely — that the broad 
velocity component is tracing outflowing material which, combined 
with extinction, gives rise to the consistent 50 km s“^ blueshift 
with respect to the disk component. But our measurements sug¬ 
gest that the outflowing gas clumps traced by H a in NGC 7552 
have not yet been entrained in the larger wind. They may be at the 
base of the wind, not yet having been accelerated to typical out¬ 
flow velocities, or may be cluster-scale outflows still embedded in 
the disk. In either case, it appears difficult to recover the true bulk 
flow rate from the observed H a emission profile in the presence 
of turbulent motions and extinction. Furthermore, the kinematics 
of the H a profile are not representative of the high velocities that 
gas entrained in the galactic wind eventually reach. As a result, UV 
absorption lines seem to more-completely trace the warm ionized 
galactic wind of NGC 7552. 


7 CONCLUSION 

We have analyzed spatially-resolved, medium-resolution, optical 
spectra of the barred spiral galaxy NGC 7552, using the emission 
lines [N II] AA6548,6583; H a; and [S II] AA6716,6731. We have 
also analyzed HST/COS far-UV spectra along a sight-line coinci¬ 
dent with our visible observations. We summarize our main results 
as follows: 

(i) We detect two velocity components in the optical spectra. We 
measure a narrow emission component (FWHM ~ 100 km s“^) 
that is consistent with H II regions in other LIRGs. We detect an 
additional broad component that is typically blueshifted by | At)| = 
30 - 50 km s“^and shows line widths of 240 - 350 km s“^ FWHM. 
Both components exhibit high electron densities (rie > 500 cm“®). 

(ii) Using methodology similar to other studies, we use the 
broad emission component to measure a mass outflow rate in H a 
of Mout = 5.0lg g Mq yr“^. This implies a mass loading factor 
of 77 = 0.3 - 0.5. These values have been corrected for an extinction 
of E{B — V) = 0.64 measured from UV stellar continuum fitting 
and assuming a |Calzetti et al.|p000^ extinction law. We caution, 
however, that we do not believe the kinematics of the H a emission 
in NGC 7552 are consistent with large-scale mass outflow. A better 
estimate of the outflow velocity yields Mout = 2.2 Mq yr“^ (77 = 
0 . 1 - 0 . 2 ). 
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(iii) Line ratio diagnostics suggest that the H Q-emitting gas 
could he at least partially shock-ionized (log([S II]/H q) « —0.5). 
Such shocks may be driving the large line widths of the broad ve¬ 
locity component. Turbulent broadening may also contribute to the 
large line widths. These additional broadening terms suggest that 
the line width is not representative of the outflow velocity (i.e., 

Wout 7^ Aw max)- 

(iv) We independently measure the mass outflow rate using UV 
silicon absorption lines and find a rate of M„ind = 19 M© yr“^ 
and an implied mass loading factor of ?7 = 1.3 - 1.9. This value 
is a factor of 8.5 higher than our best estimate in H a. Our mea¬ 
surements may underestimate the gas column density N(H), which 
would lead to an underestimate of the UV-measured Mwind by fac¬ 
tors of a few. 

(v) Although the Ly a profile is consistent with a 50 km s“^ 
cluster-scale spherical outflow, such an outflow does not account 
for a large amount of high-velocity absorption. This high-velocity 
absorption is consistent with being material entrained in a galactic 
wind. 

Given the vigorous star formation on-going in NGC 7552 
(ATsfr ^ 15 Mq yr~^), the galaxy is likely producing a 

supernova-driven galactic wind. It is not clear, however, that the 
H a emission is an adequate tracer of this wind. The optical emis¬ 
sion lines show evidence of being shocked, suggesting that the 
broad linewidths are not the result of bulk gas flows, and out¬ 
flow velocity estimates such as Awmax,2 may be substantial over¬ 
estimates, leading to over-estimates of the mass outflow rate. The 
high electron densities of this broad component also suggest it ex¬ 
ists at small radius/disk height (< 1 kpc), rather than at ~ 5 kpc as 
suggested by some authors. A large column of UV-absorbing gas 
exists at high outflow velocities (Aw > 400 km s“^) that is un¬ 
detected in our visible spectra, suggesting that the H a is tracing 
material that is either in the inner regions of the galactic wind or is 
tracing local cluster outflows that may never escape the disk. The 
Ly a emission profile appears to trace similar regions, but strong 
absorption at high blueshifted velocities of both Ly a and of low- 
ionization metal species shows that UV absorption is a more sensi¬ 
tive tracer of galactic outflows. Accurately measuring the H a out¬ 
flow mass using this methodology likely requires substantial cor¬ 
rections to account for this undetected high-velocity gas. 
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